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Abstract 

We describe a solving semi-decision method based on examination of 
the rational structures of the generalized integrating factors of first-order 
ODEs. We propose a conjecture that for some family of equations of the 
type dy/dx = P{x,y)/Q{x,y), with P and Q polynomials only in y (or in 
x), the general form of the structures of generalized integrating factors are 
rational in y (or in x). In such a way one can obtain a differential-algebraic 
polynomial system for undetermined parameters of the structures. The 
successful solution of this system (it is sufficient to find any particular 
solution) automatically leads to Ending the general solutions of ODEs. 


1 Introduction 


Equations of the type 


dyjx) 

dx 




( 1 ) 


where f{x,y) = P{x,y)/Q{x,y) with P and Q polynomials in y arise in many 
different physical fields. There are some methods for finding the solutions for 
such type of equations. First of all it is the famous method of integrating factor. 
A remarkable method based on the knowledge of the general structure of the 
integrating factor was developed by Prelle and Singer ||] and Singer Q for the 
case when P and Q are polynomials in both arguments x and y. 

There are extensions of this method in some directions including ap¬ 

proaches to solve ODEs with some transcendental or algebraic terms (heuristi- 
cally)) @-§. 

Roughly speaking, the Prelle-Singer type methods consist of three parts. On 
first of them it is determined (or supposed) the general structure of the inte¬ 
grating factors for some family of ODEs or for some type (elementary, algebraic, 
Liouvillian and so on) of solutions. On the second part for given rational ODE 
one have to choose the more specific structure of integrating factor from its gen¬ 
eral form (by specifying degree bounds for polynomials in rational, exponential 
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and so on parts of integrating factor). And then on the third part one have to 
try to find the undetermined parameters of conjectural structure. 

The distinctive feature of the Prelle-Singer type of methods is that (except 
some relatively simple cases) we are not able theoretically define a degree bounds 
for polynomials of the structures of integrating factors, so we can not miss 
vexatious second step. 

In present paper we propose an approach based on examination of structures 
of generalized integrating factors of equations with P and Q polynomials only 
in y (or only in x). 

The paper is organized as follows. In section 2 we introduce the generalized 
integrating factors. We note in Section 3 that rational hypothesis describes not 
only rational equations. Section 4 is devoted to consideration of general struc¬ 
tures of the generalized integrating factors for rational ODEs. Here we introduce 
’’Order guides” which is helpful tool for defining of procedure parameters. The 
important moment of the method - simplification of the system for unknown 
coefficients is discussed in Section 5. In Section 6 we give formulae to express 
the equation solutions through the generalized integrating factor. One of the 
ways to generalize structures for rational ODEs is considered in Section 7. In 
sections 8 and 9 we briefly touch on the cases when P and Q polynomials in x 
only, and in both x and y. 


2 Generalized integrating factors 


In we demonstrated on examples the solving abilities of the method based 
on structures of solutions of linear first-order PDE associated with (|^) 


dC{x,y) 

dx 


f{x,y) 


dC{x,y) 

dy 


= 0 . 


We examined there the structures of the type 


( 2 ) 




Po{x,y) Pt{x,y) 

qo{x,y) ''®(^,?/)’ 


( 3 ) 


or equivalently 


C{x,y) = exp( 


Po{x,y) 

qo{x,y) 


-.yr ( Pi{x,y) 


r, 


( 4 ) 


where po, qo, Pi, qi are polynomials in j/; ai are constants. The structures of type 
(i, i) are complicated enough. They may have many terms and we have to 
look over many preliminary propositions on the second step in choosing orders 
of polynomials po, qo, pi, qi and a^. 

But from (i we can note that this structure is the result of integration of 
a rational function in y. It brings us an idea that we can find such character¬ 
istic, which lead to easy but more general structures and which has solvable 
connection at least with ^(x,y). 
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It is easy to see that for classical integrating factor My = , which is deter¬ 

mined by the following equation 


^ dffiy ^ 
dx dy 


(5) 


the possible structures have the same form as (||) although here they circum¬ 
scribe a more rich ODEs family. 

The same structures has the second integrating factor y,x = ^ with the 
following equation 


d My 
dx f 


djh 

dy 


= 0 . 


( 6 ) 


Let us now introduce new characteristic 


_ d dC 
dy^^ dy 


(7) 


We will call such type of characteristics as generalized integrating factors. 
The Myy is determined by the following non-homogeneous equation 


df^yy dfnyy ^ 

dx dy dy"^ 


( 8 ) 


Here we can see that for some rational f{x,y) in y there exist a particular 
solution for fXyy in form of rational function in y. 

If we know any particular solution for Myy we can find both integrating factor 
My and PDE solution ^ and finally the general solution of equation i) y(x) (see 
Section 6). 

It is obvious, that rational conjecture for fiyy gives as particular cases all 
forms, considered above, for My and f. 

The same observation may be done for the following generalized integrating 
factor 


l-^yx — 


-^In — 
dy dx' 


(9) 


The Myx is determined by the following equation 


djiyx dfy,yx dlnf dlnf _ 1 d'^f 

dx dy dx dy f dxdy 

And also for 

d dC 

The Mxy is determined by the following equation 


( 10 ) 

( 11 ) 


d Mxy 
dx f 

And at last for 


dy^xy _ d\nf d\nf 
dy dx dy 


f-^xx — 


— In — 
dx dx 


i^ = 0. 

/ dxdy 


( 12 ) 

(13) 
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(14) 


The fjLxx is determined by the following equation 

•9 ^xx d^xx I ^ ^ 

dx f dy dx"^ f 

Each of the generalized integrating factors can be expressed through another 
one and / just as yix = —fy-y For example, 


ain/ 

Myx — h-yy + , 

(15) 

l^xy = —ffXyx 

(16) 

^ , dlnf ^dlnf 

llxx- Il^yy+ J ■ 

(17) 


Then for / rational if one of generalized integrating factor is rational the rest 
can be choosen rational too. 

We can conclude that: 

Proposition: For some family of first-order ODEs with rational f{x,y) in 
y with coefficients in the complex field C(x) there exists a particular solution 
of d) for pLyy in form of rational function in y with coefficients in C{x). The 
same is valid for pyx, pLxy and fixx- 

3 Some remarks on ODE families determined by 
rational generalized integrating factor 

If we solve, for example, equation (^) with respect to /, we obtain that 

f{x,y) = [Fl{x) + J {F2{x) - J ^^dy) exp( J yyydy)dy] exp{- J Pyydy), 

(18) 

where ^ 1 ( 0 :) and F2{x) are arbitrary functions. 

It is easy to see here that for arbitrary rational pLyy f may not be an rational 
function. Even in simplest case when pyy = const f is not rational. 

When we consider structures for f in form (^) we always receive / rational 
but not all rational / are satisfied by such a form. Analogous structures for piy 
or pLx no longer give only rational /’s but they are not complete all rational / 
too. We hope that rational generalized integrating factor can complete more 
rational /. 

From another side nothing prevent us to try to find solutions of equations (^, 
(0) and ( p^ with non-rational / starting from hypothesis that respective 
generalized integrating factor is rational. 

It is interesting to note, for example, that the ’’degenerate” case when pLxy = 
0 ( or pLyx = 0) corresponds to ’’separable” ODEs. 

In the sequel however we reduce ourselves to consideration of rational /’s 
only. 
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4 Order guides for rational generalized integrat¬ 
ing factors 


Let now f{x,y) = P{x,y)/Q{x,y) and fiyy{x,y) = X{x,y)/Y{x,y), where P, 
Q, X and Y are polynomials in y. Here we have to consider P and Q as specified 


in advance, and X and Y as unknowns. That is we are going to seek fj,yy in the 


form of 

^ Xix,y) ^ T,fJoO'Uix)y^ 
Y{x,y) J:lla2,ix)y^ 



(19) 


with unknowns ali(a;) and a2i{x), N1 and N2. Then from equation we 
obtain that X and Y have to satisfy the following equation (partial derivatives 
are denoted by standard subscripts) 


Y\PQQyy - 2PQl - PyyQ^ + 2PyQQy) + Q^iXYa, - X,Y) + 

+ PQ^XYy - XyY) + QXY{PQy - PyQ) = 0. (20) 


It is a polynomial equation of the type 

N 

^H,(x)y*=0, (21) 

i=0 


where Ai{x) depend on polynomials coefficients and their derivatives. This 
equation disintegrates on a polynomial differential-algebraic system of IV -|- 1 
equations for polynomial coefficients {Ai{x) = 0}, which determines the ’’auxil¬ 
iary” system for unknown coefficients of X and Y. To solve this system we first 
of all have to capture orders of unknown polynomials X and Y. 

Let us denote orders of polynomials in y as Np^, ^Xy, ^Yy We see 
that set of orders of polynomials in equation (p^ ) is as follows 


{Np^ +2 Nq^ +2NYy -2; +Nxy P^Yy ; Np^ +^Xy +NYy-l}. (22) 

Examination of equation (|^ shows that for substantial solutions orders at 
least of two polynomials there must be equal. It leads to systems of following 
type 


Np^ + 2iVQj, + ‘2NYy - 2 = Np^ + 2 Nq^ + Nxy + Nxy - 1; 

Np^ + ‘^^Qy + ^Xy + Nxy - 1 > + Nxy + Nxy ■ ( 23 ) 

The solution of this system is Np^ — Nq^ > 1 ; Nxy — Nxy = 1 . 

The number of equations in system Nsys here is Np + 2 Nq + 2Ny — 1 = 
Np^+2Nq^ +Nxy P^Yy and number of independent unknown coefficients Njjnc 
(taking into account possible cancellation of rational function on one nonzero 
coefficient) is Nxy -I-1Vyj^-|- 1. So we conclude that Nsys — Nunc = Np^+2Nq^ — 1 
and generally speaking (except for Np^ = 0,1, Nq^ = 0 ) Nsys — Xjjnc > 0, i.e., 
there is shortage of unknown coefficients. So a solution of the system exists only 
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if some conditions (integrability conditions or constrains) for known coefficients 
of P and Q are hold true. 

Extremely appreciably here that the estimation of maximal number of con¬ 
strains maxNcons = Nsys — Nunc is not depend on Nxy and Nyy at all. Though 
contents of the constrains strongly depend on these parameters. 

We conclude that as in original version of the Prelle-Singer method we are 
not able here on the base of orders of polynomials define their common degree 
bound. For many ODEs the degree bound is completely defined by interference 
properties of all coefficients of P and Q. 

But it is obvious that the shortage of unknown coefficients here signals first of 
all that rational generalized integrating factors do not complete the whole family 
of rational first-order ODEs. 

Considering all possible systems for orders similar to (|2^) we can arrange 
the table of cases for rational piyy (Table 1.). 


Case 

NPy - NQy > 1 

Npy - NQy < 1 

1 

<o 

II 

Conditions on 

II 

a 

1 

II 

1 

NVy -NXy<l 

Nvy and Nxy 


= NQy-Npy+2 


Nsys 

Npy 2NQy-\- 

^NQy + 

Npy + 2NQy + 


+ NXy + NVy 

+ NXy +NYy+l 

+ NXy + NVy 

maxNcons 

Npy p 2NQy - 1 

3NQy 

3Npy - 3 = 3NQy 


Table 1. Order guide for rational pyy. 

Although mathematically all rational generalized integrating factors are equiv¬ 
alent, this does not mean that the actual theirs computations are equally easy. 

Similarly we can arrange the tables of cases for rational pLxy and Pxx, 
(Tables 2-4.). 


Case 

Npy - 

>1 

Npy- 

-NQy<l 

Npy - 

-NQy = l 

Conditions on 

NVy- 

-NXy = 

NVy- 

-NXy = l 

Nvy- 

-NXy<l 

Nvy and Nxy 

= Np.^ 

y -^Qy 





Nsys 

mpy 

+ NQy + 

^NPy 

+ ‘^NQy + 

‘^NPy 

+ 2iVQ„ + 


+Nx, 

y+NVy 

+ NXy 

+ NYy + l 

+Nx, 

.+NYy+l 

maxNcons 

mpy + NQy-1 

2Np^ 

+ 2No 

1 Wy 

^Npy- 

2 = 4fVQ„ + 2 


Table 2. Order guide for rational p,yx. 
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Case 

Np^ 

-Nq^>1 

NP^- 

-Nq^<1 

Np^ 

=1 

Conditions on 


1 

It 

o 


-Nx,= 

Ny, 

- < 0 

Ny^ and Nx^ 



= ^Qy 

,-Np^+l 



Xsys 

2Np, 

+ 2No + 

Np^ 

+ 3iVQ„ + 

2Np^ 

, + 2A^q„ + 


YNx, + Ny^ 


+ + 1 

+Nx^ + Ny^ 

maxNcons 

to 

+ 2iVQ„ - 1 

Np^ 

+ 3iVQ„ 

4iVp„- 

3 = 4iVQ„ + 1 


Table 3. Order guide for rational 


Case 

NPy - > 1 

NPy - ^Qy < 0 

Np^ 

-NQy=l 

Conditions on 

Ny,-Nx = 

Ny, -Nx,=0 

NYy 

-Nx^<0 

Ny^ and Nx^ 

t-H 

1 

1 

II 




Xsys 

3iVp„ + 

2Npy + Nq^Y 


3NpY 


+Nx^ + Ny^ 

YNxy YNy^Yl 

YNx^ Y Ny^ 

maxNcons 

mp^ - 1 

2iVp„+iVQ„ 

3Np^- 

■ 1 = 3iVQ„ + 2 


Table 4. Order guide for rational fj.xx- 


For some cases there are possibilities to guess a part of structure of rational 
generalized integrating factors. As a result it reduces an order of polynomial in 
(il) (or, rather, the orders Nx^ and Ny^) and somewhat simplifies the problem. 
Let us consider as an example fjLyy{x, y) = X{x, y)/{Q{x, y)Y(x, y)), where Q is 
the denominator of /. Here we get an analog of equation (|^) in the following 
form 

Y^[PQQyy - 2PQl - PyyQ^ + 2PyQQy) - Q^iXYx + Xx)+ 

+ PQiXYy - Xy) + X{2PQy Y QQx - PyQ) = 0 (24) 

with the following table (Table 5) 
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Case 

NPy- 


NPy-NQ^<l 

Np -No =1 

JTy 

Conditions on 

NVy 

-Nxy = 

NYy-NXy = 

Nyy - Nxy < 2 - Np^ 

Nyy and Nxy 

1 

- ^Qy 

1 

II 


Nsys 

NPy 

+ NQy + 

^NQy + 

Npy +Nq^ + 


+Nxy + Nyy 

+Nxy +Ny^+l 

+Nxy + Ny^ 

maxNcons 

Np - 

y 

^NQy-l 

‘^NQy 

2Np^-2 = 2Nq^ 


Table 5. Order guide for rational ^yy = X/{QY). 


We can see from Table 5 that here the estimation of maximal number of 
constrains maxNcons = Nsys — Njjnc results in less value. 


5 Simplification of the system for unknown co¬ 
efficients 


The most natural way for finding unknown coefficients of the structures of ra¬ 
tional generalized integrating factors is solving of the system for unknown co¬ 
efficients, which follows from polynomial equations of type For example, 

for ^yy this system is determined by the equation (20). 

Of course, as in every variant of Prelle-Singer procedure from the beginning 
we have to assign a bound to the Nxy (or iVy^). The order of Nyy (corre¬ 
spondingly Nxy) is defined by ’’Conditions on Nxy and Nyy'" from the Tables 
1-5. 


As we stated above such systems are polynomial ODE-algebraic systems. On 
the first step here we have to simplify this system by any of existent differential 
elimination procedures 0, 00- The systems described here are not easy by 
no means. And it is crucial moment of the method. Even for Nsys = 5...6 the 
procedure sometimes may hang-up for days. 

If we consider ODE with symbolic parameters then under properly (by ” Con¬ 
ditions on Nxy and Ny^ ”) assigned Nxy and Ny^ and when Nsys is sufficiently 
small we always receive a simplified system with some (or without) conditions 
(constrains) for ’’known” coefficients of P and Q. Then here the ODE solution 
can be found if these constrains are satisfied and we are able to solve simplified 
system. 

In case when ODE is specified without symbolic parameters then we can 
receive or ’’System is inconsistent” (and we may augment the bound to the 
orders of unknown polynomials) or simplified system without any constrains 


have to note here that for the problems described in this paper the Rif procedures in 
Maple lead to goal faster. 











(here we are disquieted if the solution of simplified system can be expressed in 
some way). 

Despite the foregoing and the fact that it is as before a semi-decision proce¬ 
dure, the method described here can solve plenty of unsolved equations. 

Compared with pattern matching techniques the method proposed here will 
certaily resource consuming and so not be so efficient in some cases, but it 
provides a systematic way of finding solutions without a prior knowledge of 
what type the equation is. And what is more, in many cases we can use (to pick 
up speed of the process) pattern matching techniques too with precomputed (by 
described above prosedures) lists of integrability conditions. 

We are not going to give examples here. They take up a lot of space but 
one easily can generate them, e.g., by means of short Maple program like this 
example for simplest nontrivial ODE: dy/dx = l/{Yl{x)y -I- yO(a;)) 
restart; 

with (PDEtools): 

N :=0: 

n[l] ■= N : 

n[2] ■=N + 3: 

f := l/{Yl{x) *y+ Y0{x)) : 

mu := sum{a\[i]{x) *y^,i = 0..n[l])/sum{a2[i]{x) *y'^,i = 0..n[2]) : 
sys := {coeffs{collect{numer{factor{diff{f, ‘$‘(y, 2)) -|- diff{mu, x)+ 
+diff{mu, y)* f + diff{f, y) * mu)), y), y,' t')} : 
furies := {seq{al[i]{x), i = 0..n[l]), seq{a2[i]{x), i = 0..n[2])} : 
casesplit{sys, funcs, rif); 


6 The solution of equation assembling 

If we have found a solution of the system for unknown coefficients (that is we 
have found, for example, the y,yy) then it remains only to assemble the solution 
of equation (||). 

We have from (|^) that 


Q{x,y) = Fx{x)+F 2 {x) 


gj Uvvi^,v)<iy 


(25) 


where ^ 1 ( 2 :) and F 2 {x) have to be determined. If we substitute ( p5|) to (||), we 
can find that 


F 2 {x) = C 2 exp{- J Hyy dy} exp{- J (^ + / a 


,j)dx} 


(26) 
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and 


Fi{x) = C 1 -C 2 J exp{- J fiyy dy} ex-p{- J + fyyy)dx}x 

^ I I 


^</ 


^dy 

OX 


(9y 


+ fd-yy} + f exp{ J iJ-yy dy}} dx, (27) 


where Ci and C 2 are arbitrary constants. Here we may impose that Ci = 0 and 
C 2 = 1 . 

It is easy to verify that Fi {x) and F 2 (x) are really functions only in x if y,yy 
is a solution of (§). 

So, all components of C in are known now and we can obtain the general 
solution of equation (|^) in implicit form as follows 


C{x,y) = C, 


(28) 


where C is an arbitrary constant. 


7 More general structures 


We have concluded above that rational generalized integrating factors do not 
fit all rational hrst-order ODE. 

There are some possible ways for extension of forms of generalized integrating 
factors. The most direct of them based on the following idea. 

The solution of linear first-order non-homogeneous PDE, e.g., 


ox oy 


(29) 


(where / = /(x, y), g = g(x^ y)) can be expressed through particular solution of 
homogeneous PDE 


dm f dgo 
dx ■' dy 


= 0 


(30) 


by the following expression 


g{x,y) = F[r]o{x,y)] 



g{^, RootOf[r]o{^, Z) - ? 7 o(a;, y)\)d^, 


(31) 


where ^[ 5 ] is an arbitrary function, RootOf[(l){Z)\ is a root of algebraic equation 
diZ) = 0 with respect to Z. 

Further we can suppose (on the same reasons as explained above) that 
^ In is rational, then for 5 = (and F[s] = s) from (^ and ( ^l| ) we 
obtain a structure for g,y = exp{r]), which obviously is more general. 

Unfortunately this way leads to sufficiently complicated structures to be 
analyzed yet. 
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8 Rational ODEs only in x 

It is obvious that we can consider by the same way first-order ODEs with / 
rational only in x. The corresponding order guides are mirrors of Tables 1-4 
under replacement oi P Q and x ^ y. 


9 Rational ODEs in both X and y 


It is a classical case of the Prelle-Singer method. We can combine different order 
cases and arrange the order guides too. 

Here we seek solution for generalized integrating factors in form 


l^yy ~ 




E iVjS 

i=0 




i i 


(32) 


where Oy and 6^ are constants. 

Therefore we obtain the ’’auxiliary” system in form of algebraic polynomial 
system. The maximal number of constraints here depends on Nx ^, Nxy, Ny^, 
and Nxy but we are not able define a common degree bound for polynomial 
orders too. 

10 Conclusions 

We presented the semi-decision method for solving rational first-order ODEs 
based on hypothesis that generalized integrating factor of the given ODE is 
rational function. 

Besides first-order ODE the method is suited after some modification for 
solving some non-linear equations of higher order too. 
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